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The development of innovative terahertz (THz) imaging systems has recently moved in the focus of 
scientific efforts due to the ability to screen substances through textiles or plastics. The invention 
of THz imaging systems with high spatial resolution is of increasing interest for applications in the 
realms of quality control, spectroscopy in dusty environment and security inspections. One of the 
main restrictions of current THz imaging systems is the low spatial resolution which is limited by a 
lack of THz lenses with strong focusing capabilities. Here we present the design, fabrication and the 
measurement of the optical properties of spectrally broadband metamaterial-based gradient index 
(GRIN) lenses that allow one to focus THz radiation to a spot diameter smaller than the wavelength. 
Due to the subwavelength thickness and the high focusing strength the presented GRIN lenses are 
an important step towards compact THz imaging systems with strongly improved spatial resolution. 



THz physics and metamatcrials have raised a great 
deal of scientific interest in the last few years. While 
recent advances in the THz technology were strongly 
driven by efforts to transfer the fundamental knowl- 
edge about THz radiation to industrial applications, 
the metamaterial science - as the younger of both 
research fields - was mainly dedicated to the inves- 
tigation of the fundamental physical properties of 
metamatcrials. This especially applies to metama- 
terials at higher frequencies where high absorptive 
loss [1] and the enormous demands on sophisticated 
three-dimensional nano-fabrication techniques [2, 3] 
obstructs the development of competitive and appli- 
cable metamaterial-based optics. In respect thereof, 
the limitations lead to the pursuit of many interest- 
ing basic research topics such as e. g. metamatcrials 
with gain [4-7] or the investigation of quantum effects 
in metamaterials [8-10]. In contrast, metamaterials 
at microwave frequencies are advantageous in terms 
of low optical loss and ease of fabrication. In con- 
sequence, microwave metamaterials have been devel- 
oped to a maturity level that enables the integration 
into applied systems [11]. 

Another promising frequency band for the applica- 
tion of metamaterials is the range from 0.1 to 10 THz, 
the so-called THz range. Whereas THz physics was 
mostly concentrated on fundamental research about 
7 years ago, the rapid development of the technol- 
ogy and the fabrication of improved THz emitters 
and detectors lead to the opening up of new markets 
and application fields in the last few years. The spe- 
cific property of THz radiation to penetrate through 
most dielectrics as e. g. textiles, paper, concrete etc. 
has been exploited in many scenarios that are related 
to non-invasive quality control through packages [12- 
14], chemometrics of pharmaceutical substances [15], 
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terahertz imaging systems [16, 17] and security in- 
spections [18-21]. The rapid expansion of industrial 
target applications is attended by increasing demands 
for high quality optical components for the THz tech- 
nology. Moreover, many potential applications can- 
not yet be addressed by the THz technology since 
the required components are not available. However, 
such optics cannot be easily devised since the insuf- 
ficient electromagnetic response of most dielectrics 
prevents the development of standard optics that are 
based on the same principles as conventional optics 
used in laser systems. In this context, metamateri- 
als have already proven to offer a novel approach for 
modulating THz radiation [22-25], perfect absorbers 
with dcsignablc properties [26, 27], frequency filters 
[28] or wave plates [29]. Another application field 
where metamaterials can provide an important con- 
tribution is the THz imaging technology. One of the 
main restrictions of current compact state-of-the-art 
THz imaging systems is the low spatial resolution 
which is limited to about 1 mm due to a lack of avail- 
able THz lenses with strong focusing capabilities. It 
is obvious that the optical resolution of such imaging 
systems could be significantly improved if lenses with 
strong focusing capabilities and low aberrations were 
available. 

Here we present the design, fabrication and the 
experimental testing of a three-layer metamaterial- 
based GRIN lens that allows one to focus THz radia- 
tion to a spot diameter smaller than the wavelength. 
We also compared the optical properties of the 3-layer 
GRIN lens to the focusing behavior of a 1-layer GRIN 
lens that is significantly easier to fabricate. We could 
show that a 1-layer lens already focuses THz radia- 
tion to a spot size in the order of the wavelength of 
the THz beam. The lenses offer broadband opera- 
tion, are very thin and allow aberration-free optical 
imaging due to the avoidance of curved boundaries 
in the beam path. In contrast to most other materi- 
als, metamaterial lenses possess the potential to serve 
as adaptive optical components since they allow one 
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Figure 1. (a) Structural design of the unit cell. The underlying geometry was based on annular slots in a copper plane, 
(b) The refractive index of the metamaterial was changed by varying the inner radius of the annular slot between r = 18 
and 23 urn. This resulted in a total change of the refractive index from 0.08 to 1.65 at a frequency of 1.3 THz. (c) 
Microscope picture of the 3-layer GRIN lens, (d) Refractive index in dependence of the radial distance from the center 
of the GRIN lens. 



in principle to tune the optical properties by optical 
or electronic means. In consequence, the presented 
GRIN lenses can be considered as an important step 
towards compact THz imaging systems with strongly 
improved spatial resolution. 

Fig. 1(a) illustrates the underlying structural de- 
sign of the GRIN lens [30, 31]. The unit cell geom- 
etry was based on annular slots in a copper plane. 
Thereby, the refractive index of the unit cell could 
be strongly varied by changing the inner radius of 
the annular slot between r — 18 and 23 um at a con- 
stant slot width of 3um. The copper structure was 
fully embedded in a benzocyclobutcne background 
matrix with a permittivity of e = 2.67 and a loss 
tangent of S = 0.012. The dimensions of the unit 
cell were 60 2 x 40 pm 3 . The copper layer had a 
thickness of 200 nm and an electric conductivity of 
cr = 5.8 x 10 7 S/m. We retrieved the refractive in- 
dex of the specific unit cell designs by use of three- 
dimensional full wave simulations with the software 
CST Microwave Studio. The numerical calculations 
showed that the investigated structure supports a re- 
fractive index gradient between the minimum and 
maximum values from 0.08 to 1.65 at a center fre- 
quency of 1.3 THz (see Fig. 1(b)). For the deter- 
mination of the extrema we restricted ourselves to 



the effective medium regime and took carefully the 
boundary to the Bragg regime into account since op- 
eration of the lens in the Bragg regime would lead 
to undesired spatial patterning of the focus and thus 
distort the optical quality of the focusing device. For 
the lens design we chose a radially symmetric refrac- 
tive index gradient. For this purpose, we introduced 
a spatial variation of the refractive index by arraying 
unit cells of different inner slot radius such that the 
refractive index gradually decreased from the center 
of the GRIN lens. This is indicated in the microscope 
picture in Fig. 1(c) and Fig. 1(d). In more detail, the 
spatial dependence of the refractive index could be 
approximated by a parabolic index profile. Such a 
lens is conceived to focus incoming radiation since it 
is expected from the refractive index profile that an 
electromagnetic wave travels faster in the outer an- 
nular segments of the lens than in the inner annular 
segments. In this regard the refraction power corre- 
lates with the maximal achieved difference between 
the phase advance of the wave across the outer bor- 
der and through the center of the GRIN lens. Thus, a 
thicker lens should result in a higher refraction power. 
For this reason, we compared the optical performance 
of a 1-layer and a 3-layer thick GRIN lens. 

We fabricated both the 1- and the 3-layer GRIN 
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Figure 2. Schematic of the experimental setup for the beam profile measurement and the determination of the focus 
position of the GRIN lens. The THz electric field behind the lens is detected by exploiting the electro-optic effect 
in gallium phosphide (GaP). The optical probe beam is scanned over the THz beam to determine the spatial field 
distribution with a spatial resolution of 60 pm that is prescribed by the spot size of the optical beam. PM1,2: parabolic 
mirrors, LI: lens for the optical beam, Ml, 2: dielectric mirrors for the optical beam, A/4: quarter-wave plate 



lens by UV lithography. That way we obtained a 
1-layer GRIN lens with a thickness of 40 um and a 
3-layer lens of 120 um thickness. We investigated the 
GRIN lenses with respect to their optical properties 
by mapping the spatial distribution of the complex 
electric THz fields. A schematic of the used experi- 
mental setup is shown in Fig. 2. We generated short 
THz pulses with a duration of 2 ps and a spectral 
bandwidth from 0.1 to 2.0 THz by focusing ultra- 
short laser pulses with a duration of 15 fs between the 
poles of a photoconductive switch. The emitted THz 
beam was polarized along the x-direction. In order 
to match the THz beam size to the aperture of the 
investigated GRIN lenses, we collimatcd the emitted 
THz pulses by an off-axis parabolic mirror PM1 and 
focused the THz beam to a spot diameter of about 
1.5 mm by a second parabolic mirror PM2. Hereby 
and in the following the spot diameter is defined by 
the full width at half maximum (FWHM) of the in- 
tensity distribution. We carefully aligned the GRIN 
lens with respect to the optical axis of the THz beam 
and positioned it in the beam waist of the THz ra- 
diation. By this method, we assured that the wave 
fronts of the incident THz radiation were plane at 
the entrance facet of the GRIN lens. Uttermost care 
was advised for the lateral alignment of the GRIN 
lens since the THz waves would have experienced 
an asymmetric refractive index gradient if slight lat- 
eral shifts of the center point of the GRIN lens had 
been introduced. Such a misalignment and the result- 
ing asymmetric refractive index profile would cause 
highly undesirable beam shape distortions. 

To determine the beam profile and the spot size 
of the focused THz beam we measured the spatial 
distribution of the complex electric field by use of 



electro-optic sampling in a reflection geometry [32, 
33]. This configuration is illustrated in Fig. 2 and 
the measurement technique is described in detail in 
the methods section of this letter. 

In order to identify the position of the focal plane 
and the focus diameter of the 1- and 3-layer GRIN 
lens we measured the transversal x-y-intensity pro- 
file of the THz beam in dependence of the position 
z along the propagation direction. Examples of the 
obtained 2D intensity profiles are shown as insets of 
Figs. 3(a) and 3(b). From the 2D profiles we deter- 
mined the ID intensity distribution of the THz beam 
along the cross section line in y-direction through the 
center of the intensity distribution. This is illustrated 
in Fig. 3(a) for the 1-layer lens and in Fig. 3(b) for 
the 3-layer lens as measured at a center frequency of 
1.3 THz. Thereby, the position z was defined relative 
to the position of the corresponding focal plane. From 
the ID intensity distributions we extracted the spot 
diameter D in dependence of the z-position as illus- 
trated in Fig. 3(c) and Fig. 3(d) for the 1- and 3-layer 
lens, respectively. For comparison, we also plotted 
the spot diameter that we obtained from numerical 
full wave calculations. As can be seen from Figs. 3(c) 
and 3(d), significant deviations of the theoretical spot 
diameter calculations from the experimental data are 
observed for the 1-layer lens whereas the agreement 
between the experimental and theoretical curves is 
very reasonable for the 3-layer lens. This was ex- 
pected, since the 1-layer lens is more prone to surface 
effects that vary strongly with residual imperfections 
of the fabricated samples. In contrast, the 3-layer 
lens revealed bulk properties and the impact of sur- 
face effects on the electromagnetic properties of the 
lens was only minor. 
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Figure 3. (a) ID intensity profile of the THz beam for different z-positions in propagation direction of the THz wave 
for the 1-layer GRIN lens. The intensity profiles were obtained by extracting the intensity values along the ID cross 
section line of the 2D transversal x-y-intensity profiles along the y-direction as shown in the inset. The z-position was 
measured relative to the focal plane, (b) Same as (a) for a 3-layer GRIN lens, (c) THz beam diameter D as defined 
by the full width at half maximum (FWHM) of the intensity profile of the beam normalized to the wavelength A in 
dependence of the z-position in propagation direction of the THz wave for the 1-layer GRIN lens. For comparison, the 
theoretical data as obtained from 3D full wave simulations are plotted as solid lines, (d) Same as (c) for a 3-layer GRIN 
lens. 



As expected, both GRIN lens designs provoked a 
strong focusing of the THz radiation. By interpola- 
tion, we determined the focal plane of the 1-layer lens 
to be located 1.1 mm behind the exit facet of the lens. 
From the measurements we determined a focal beam 
diameter of D — 320 um at 1.3 THz which relates to 
D/X = 1.45 where A is the wavelength of the THz 
radiation. For the 3-layer lens, the focal plane was 
located 0.8 mm behind the lens exit surface. We mea- 
sured a focus diameter of 220 um at 1.3 THz which 
corresponds to a ratio D/X = 0.96 and thus indicates 
slight subwavelength focusing. A focus of subwave- 
length size was expected from numerical calculations 
for the 3-layer GRIN lens. The stronger focusing ca- 
pabilities of the 3-layer lens in comparison with the 
1-layer lens can be intuitively understood by consider- 
ing the focus mechanism of a GRIN lens. The longer 
the propagation length through the GRIN lens is, the 
more the phase difference between the outer part and 
the inner part of the propagating wave accumulates. 
For this reason, we expect a larger curvature of the 
phase fronts of the beam which implicates a stronger 
focusing. 

We also studied the spectral operation bandwidth 
of the lenses by measuring the diameter D of the THz 
beam normalized to the wavelength A at the focal 



plane for different center frequencies (see Fig. 4). For 
comparison, we plotted the numerically calculated ra- 
tio D/X of the considered lenses in dependence of 
the center frequency. Equivalent to the foregoing 
discussion, the numerical and experimental results 
show better quantitative agreement for the 3-layer 
lens (Fig. 4(b)) than for the 1-layer lens (Fig. 4(a)). 
The measurements evidence that the GRIN lenses 
operated in a frequency range from 1.2 to 1.5 THz. 
Within this range, the ratio D/X varied between 1.45 
and 1.75 for the 1-layer lens and between 0.94 and 
1.05 for the 3-layer lens. Therefore, we could achieve 
weak subwavelength focusing of THz radiation with 
a 3-layer GRIN lens over almost the complete fre- 
quency band for which the lens was designed. Note 
that the GRIN lenses operated in the Bragg regime 
for frequencies higher than 1.5 THz (grey shaded re- 
gion in Fig. 4). Since the beam shape was distorted 
by scattering in this frequency range the FWHM def- 
inition of the beam diameter was not longer valid. 
As mentioned before, lens operation in the scattering 
regime is not recommended. 

In conclusion, we have demonstrated the design, 
fabrication and measurement of the optical proper- 
ties of a 1-layer and a 3-layer GRIN lens and com- 
pared the optical performance of both lens types. 
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Figure 4. (a) Experimental and numerical values of D / A as defined by the ratio between the diameter D of the THz beam 
and the center wavelength A of the THz radiation for the 1-layer GRIN lens. At frequencies higher than 1.5 THz the 
lens no longer operated in the effective medium regime and the optical behavior was governed by scattering. Operation 
at frequencies higher than 1.5 THz is therefore not recommended, (b) Same as (a) for the 3-layer lens. As envisioned 
by the design, the GRIN lens focused incident THz radiation to a focus diameter smaller than the wavelength. Over a 
frequency band between 1.2 and 1.5 THz the ratio D/X was close to unity. 



We showed that a 1-layer GRIN lens already offers 
strong focusing capabilities and allows one to focus 
THz radiation to a spot diameter in the order of the 
wavelength. At a frequency of 1.3 THz we obtained 
a focus diameter of D = 320 um. Even slight sub- 
wavelength focusing was expected and obtained for a 
3-layer GRIN lens. We measured a focus diameter of 
D = 220 um at a frequency of 1.3 THz which corre- 
sponds to a ratio between the focus diameter D and 
the wavelength A of D/X = 0.96. Both lens types op- 
erated over a comparably large frequency range from 
1.2 to 1.5 THz due to the non-resonant design of the 
metamaterial structure. The subwavclcngth focusing 
capabilities of the demonstrated metamaterial-based 
GRIN lenses in a broad frequency range can be con- 
sidered as an important step towards the improve- 
ment of the spatial resolution of compact, industrial 
THz imaging systems. 



METHODS 

Measurement setup 

As an electro-optic material we employed a gallium 
phosphide (GaP) crystal with a thickness of 400 pm 
and an aperture size of 20x20 mm 2 . The GaP crystal 
was cut in the (llO)-direction such that the electro- 
optic effect was only sensitive to the x-polarization of 
the THz beam. The detection scheme was conceived 
as follows: Due to the electro-optic effect, the inci- 
dent THz pulses induce an instantaneous change of 
the refractive index of the GaP crystal. This change 
is proportional to the momentary amplitude of the 
electric THz field. Since the electro-optic effect in 



GaP is anisotropic, it is possible to determine the 
amplitude of the electric THz field by measuring the 
change of the polarization state of ultrashort probe 
pulses that propagate through the crystal. For this 
purpose, we focused linearly polarized probe pulses 
to a focal spot size of 60 pm at the position of the 
GaP crystal. The wavelength of the probe pulses was 
800 nm and the pulse duration was 40 fs. The sur- 
face of the GaP crystal facing the THz beam was 
coated by a highly reflective layer for a wavelength 
of 800 nm while the opposite facet was equipped by 
an anti-reflection coating for the identical wavelength 
range. The optical coatings were designed in such a 
way that they did not affect the propagation proper- 
ties of the THz beam. Consequently, the THz pulses 
could penetrate the GaP crystal without significant 
loss. On the other hand, we herewith ensured that 
the probe pulses could enter the GaP crystal with 
strongly reduced loss and were efficiently reflected at 
the other facet of the crystal where the THz pulses 
entered. During the propagation through the crys- 
tal the probe pulses changed their polarization state 
due to the optically induced change of the refractive 
index by the momentary THz field. Because of the 
anisotropy of the electro-optic effect, the difference in 
the induced phase delay of two orthogonally polarized 
field components of the probe pulses resulted in ellip- 
tically polarized probe pulses when the THz field was 
non-zero. Hereby, the orientation of the major and 
minor axes of the polarization ellipse was sensitive 
to the polarity and thus the phase of the momentary 
THz field. To enhance the sensitivity of the measure- 
ment we analyzed the polarization of the probe pulses 
by combination of a quarter wave plate, a wollaston 
prism and a balanced detector. Since we obtained the 
probe pulses from deviating a small fraction of the en- 
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ergy of the pump pulses used for the THz generation 
by a beam splitter, the probe pulses were inherently 
coherent with the generated THz pulses. Further- 
more, the pulse length of 40 fs of the probe pulses 
was a factor of 50 smaller than the pulse length of 
the THz pulses. This allowed us to sample the tem- 
poral shape of the electric field of the THz pulses by 
delaying the probe pulses with respect to the THz 
pulses. The coherence between the THz pulses and 



the probe pulses in combination with the phase sen- 
sitivity of the used detection scheme enabled us to 
measure both the amplitude and the phase of the 
electric field. Furthermore, since the optical probe 
beam was focused to a spot size of 60 um the spatial 
resolution of the system was only limited by the spot 
size of the probe beam. Thus, we could measure the 
spatial field distribution of the THz beam with sub- 
wavelength resolution with respect to the wavelength 
of the THz radiation. 
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